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In this work we present a model and a method to study integer quantum Hall (IQH) systems. Making use of 
the Landau levels structure we divide these two dimensional systems into a set of interacting one dimensional 
gases, one for each guiding center. We show that the so-called strong field approximation, used by Kallin and 
Halperin and by MacDonald, is equivalent, in first order, to a forward scattering approximation and analyze the 
IQH systems within this approximation. Using an appropriate variation of the Landau level bosonization method 
we obtain the dispersion relations for the collective excitations and the single particle spectral functions. For the 
bulk states, these results evidence a behavior typical of non-normal strongly correlated systems, including the 
spin-charge splitting of the single particle spectral function. We discuss the origin of this behavior in the light 
of the Tomonaga-Luttinger model and the bosonization of two dimensional electron gases. 
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PACS numbers: 71.10.Pm, 71.70.Di, 73.43.-f, 73.43.Cd, 73.43.Lp 



I. INTRODUCTION 

The integer quantum Hall effect (IQHE) 1 appears in two 
dimensional electron gases (2DEG) under a strong magnetic 
field, B, with the electronic Landau level filling factor, V, 
equals to an integer number. In this case the non-interacting 
ground state presents V completely filled Landau levels and 
the excitons generated by the electron-electron interaction in- 
volve particles in different levels. At stronger magnetic fields 
the degeneracy, Afy = S/2nl 2 (where S is the area of the 2DEG 
and / = y/ch/ (eB) is the magnetic length), of a single Lan- 
dau level exceeds the number of electrons. In this case many 
non-interacting ground states can be constructed with all the 
electrons in the macroscopically degenerate lowest Landau 
level. In this case the usual perturbative theories, which lead 
to the Fermi Liquid behavior, are invalid and the effect of the 
electron-electron interaction turns out to be much more dra- 
matic, giving rise to a new state of the matter and to the frac- 
tional quantum Hall effect (FQHE) 2 ^ 3 .. 

The study of collective excitations plays an important role 
in understanding both systems. The existence of magneto- 
rotons 4 ' 5 -^ was verified through inelastic light scattering 7 and 
phonon spectroscopy 8 ' 9 . However, in the integer quantum 
Hall systems, at very low electronic densities (less than 
10~ 10 cm~ 2 ) multiple magneto-rotons which are not accessi- 
ble to the usual perturbative theories (possibly because they 
neglect the mixing of the Landau levels) were observed 10 , 
demonstrating an incomplete understanding of these excita- 
tions. Additionally, in the FQHE the coupling between the 
particle excitations with the magneto-plasmon mode close to 
the cyclotron frequency, C0 C , seems to be the key for a Hamil- 
tonian theory of these systems 11 , and to the origin of most 
part of the notable properties of their quasi-particles, the com- 
posite fermions 1213 . Here we will show that the collective 
excitations also seem to have a strong influence in the bulk 
quasi-particles properties even in the IQHE. 

In this work we will explicitly treat the integer quantum 
Hall systems, although the method we present, together with 
the Chern-Simons theory 14 , also has applications in the frac- 



tional case 15 . Assuming a system where the Coulomb en- 
ergy, e 2 /el, is much smaller than the cyclotron energy, h(O c = 
heB/(mc), the so-called strong field approximation 4 - 5 (where 
the mixing between the Landau levels is neglected) can be 
used. Under this approximation, the collective magneto- 
plasma modes were determined by the Green's function 
method 4 and time-dependent Hartree-Fock theory 5 . In this 
work, we will see that these modes can also be obtained by 
a magneto-exciton bosonization method. First, making use 
of the natural structure of the Landau levels, we divide the 
2DEG into a set of Afy one dimensional interacting channels, 
one for each guiding center occupied by V particles. We show 
that, to first order in (e 2 /el) / h(£) c , the strong field approxi- 
mation is equivalent to a forward scattering approximation in 
this model of one dimensional interacting channels. Thus, we 
will allow for the transfer of energy and momentum between 
the different channels mediated by the direct and exchange 
Coulomb interactions, conserving the net number of particles 
in each channel. We then use a variation of the Landau level 
bosonization 16 which is ideal for a treatment of the interaction 
effects in this model. 

One of the major advantages of this bosonization method 
is that, through the Mattis-Mandelstam bosonic representa- 
tion of the fermionic operato r 17 i 18 , one is be able to deter- 
mine the single particle properties directly from the collective 
modes, differently from the usual perturbative theory. An- 
other advantage is that this Mattis-Mandelstam representa- 
tion is an identity-^ which, technicaly, can be used later to 
extend the model, including the net transfer of particles be- 
tween neighbor channels according to their relative angular 
momentum. We believe that bosonization procedures similar 
to this one can be important tools to clarify the relation be- 
tween the quasi-particle properties and the collective modes 
in the FQHE. Indeed, the Hamiltonian Theory of the FQHE 11 
is very similar to a bosonization of the degrees of freedom 
related to the magneto-plasmon excitation. 

The organization of the paper is as follows. In section II we 
review the free and the interacting Hamiltonians of a two di- 
mensional electron gas in a perpendicular magnetic field writ- 
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ten in the Landau level basis and in the symmetric gauge. In 
section III we show the equivalence between the strong field 
approximation and the forward scattering approximation and 
then describe how to treat the 2DEG, including the direct and 
exchange interactions between the electrons in a bosonization 
approach. In section IV we present a comparative analysis 
of our bosonized model with the Tomonaga-Luttinger model 
and the 2D bosonization. Next we describe the strategy to deal 
with spin-full fermions in section V. The spectrum of neutral 
excitations obtained is presented in section VI. The construc- 
tion of the fermionic operator and the determination of the 
single -particle spectral functions is accomplished in section 
VII and our conclusions are presented in section VIII. 



II. FERMIONIC REPRESENTATION 

We consider a 2DEG in a strong perpendicular magnetic 
field at zero temperature with an integer number, V, of fully 
occupied Landau levels. We start from a simpler description 
of spinless electrons and generalize the model to the case of 
electrons with spins in section VII. 

The non-interacting part of the Hamiltonian of spinless 
fermions of mass m in a plane perpendicular to an uniform 
magnetic field is given by 



#0=2^ /rfW(r) 



P--A(r) 



TO 



(2.1) 



where (r) is the fermion field operator, which can be ex- 
panded in the Landau level basis and in the symmetric gauge, 
A (r) = f z x r, as 



= £(r|n,m)c„ )W . 



(2.2) 



Here \n,m) are the eigenstates of the one particle Hamiltonian, 
where n refers to the Landau level and m to the guiding center, 
withSl 



(r\n,m) = G m , n ( i-r 

Villi 2 V / 



where z = x + iy. For m> n, 

r I \ fnl (-ii 

G min {z) = W— — 
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(2.3) 



(2.4) 



with L™ " being the generalized Laguerre polynomial. The 
Hamiltonian Hq is diagonal in this Landau levels basis 



N® — 1 oo 

H = h(O c £ £(n+l/2)4 m c„, 

m=0 n=0 



(2.5) 



Given the guiding center degeneracy, Afy, the non-interacting 
ground state is defined by uniformly filling V Landau levels of 
each guiding center m : 

ty-lv-l 

i g o>= n rRjo), 

m=0 n=0 



where Afy = S/2nl 2 is the number of flux quanta crossing the 
area of the system, 5. The neutral excitations in the non- 
interacting system are electron-hole pairs with a hole in a Lan- 
dau level p < V — 1 at the guiding center m and an electron in 
a level p + n> \ — 1 at the guiding center m'. 

In the Landau levels basis the Coulomb interaction, V (r) = 
e 2 /er, is formally given by 



Hi = - £ {n\,m\;n2,m2\V\n\,m\;n 2l m 2 ) 



{n.m} 

* i ,m i ^ ?i2 ,ni2 ,ni2 ,m\ ' 



(2.6) 



where the interaction matrix element can be written in the mo- 
mentum representation as 



(n i , m i ; n2 , n?2 1 V \ n\ , m\ ; n 2 , m 2 ) 



1 



n .. , d 2 qV{q)e- 
(271) 

xG*, (lq*)G*, (la) 

XGn 2 ,n' 2 (^q*)G m2 ^(/q), 



with V (q) = 2ne 2 /eq and G defined in d2~4l . 



III. THE FORWARD SCATTERING APPROXIMATION 
AND THE BOSONIC DESCRIPTION 

A. First Order Equivalence Between the Strong Field 
Approximation and the Forward Scattering Approximation 

Inspired by the linear energy dispersion of the quantum 
number n and the degeneracy of the quantum number m, one 
can make an analogy between the quantum Hall system and 
the problem of many interacting one-dimensional electron 
gases. In fact, the elementary neutral excitations spectrum 
in the integer quantum Hall regime, consisting of electron- 
hole pairs excited over a completely filled set of Landau lev- 
els, has been calculated using the diagrammatic Green's func- 
tion approach of Kallin and Halperin 4 and the time-dependent 
Hartree-Fock approximation of MacDonald 5 . In both meth- 
ods the excited states are constructed essentially in a single 
mode approximation: (i) considering only the presence of 
a single electron-hole pair at a time, (ii) over a ground state 
that neglects any Landau level mixing. Under these approxi- 
mations, named strong field approximation, the magnetoplas- 
mon bands were calculated including different order terms in 
the interaction strength, y = (e 2 /e/) / h(ti c <C 1. Kallin and 
Halperin include only some first order terms, not all of them, 
once they do not consider the ones associated with the decay 
of one exciton into two. MacDonald includes some second 
order terms, although also excluding the first order ones ne- 
glected by Kallin and Halperin. The effect of the decay of ex- 
citons was partially considered by Cheng 19 . In Fig.([TJ, which 
was constructed according to the conditions (i) and (ii) above, 
we show that, to first order in y, the strong field approxima- 
tion is analogous to the forward scattering approximation in 
a set of one-dimensional electron gases, where the net num- 
ber of particles in each channel (here the guiding center m) 
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is conserved. Indeed, due to the conservation of the angular 
momentum, the forward scattering approximation will be al- 
ways restricted to first order effects in y, as we discuss in Sec. 
(III.C). But, as we stressed above, this is not a very restrictive 
condition, because the previous methods of calculation are not 
even able to include all these first order effects. 

We saw above that if we formally divide a two dimensional 
IQH system into m interacting one dimensional channels, the 
strong magnetic field turns the neglect of the net transfer of 
particles between the "m" channels into a reasonable approxi- 
mation. So, this system looks very convenient for the use of a 
bosonization technique, once we can skirt the treatment of the 
transfer of particles between many different one dimensional 
channels in this first approximation. The latter is known to be 
the main difficulty with the higher dimensional bosonization 
method a 20 ! 22 . However, despite the problem with the many 
channels particle transfer, there are many techniques devel- 
oped to study the effects of the net transfer of particles be- 
tween two channels in a bosonization method. This problem 
is analogous to the inclusion of the backward scattering in 
the usual bosonization of a one-dimensional interacting elec- 
tronic system, which makes the transfer of particles between 
the right- and left-moving channels 4 ^. So, in a future work, 
we intend to included the transfer between two neighbors "m" 
channels. This neighbor transfer seems to us a promising al- 
ternative expansion parameter to the interaction strength y, as 
we discuss further in the conclusion. However, in this work, 
we will restrict ourselves to the first order terms in the strong 
field approximation or, equivalently, to the forward scattering 
approximation. 



B. The Bosonic Description 

Motivated by the above discussion, we will introduce 
bosonic operators associated with electron-hole excitations 
that conserve m, as in the bosonization scheme previously de- 
veloped by two of us 16 . If we extend the Hilbert space of all 
guiding centers to negative energy values, in the spirit of the 
Luttinger model 18 ' 24 , we get that the operator 



1 00 
^ =— V c f 

n,m r- ^n+p,. 



(3.1) 



and its Hermitian conjugate obey exact bosonic relations: 



J n,m^n' .m' 



0. 



(3.2) 



This extension is a well controlled approximation as long as 
the negative energy states are kept inert. So, the system's 
Fermi energy must be much higher than the characteristic 
Coulomb interaction, vh(O c » e 2 /el, which is essentially the 
strong field approximation. 

The Hamiltonian (12.51 ) can be rewritten in terms of the pre- 
viously defined bosonic operators as^ 



H = h(o c Y, 

m=0 



n>0 1 



(3.3) 



















(a) 





(c) 



(d) 



Figure 1: (a) A single exciton created over a ground state that com- 
pletely fills the v — 1 lowest Landau levels, neglecting any Landau 
level mixing (ii). Consider now only first order processes, where 
the total energy is conserved, An; = 0. If we neglect multiexciton 
states (i), processes like (b) or (c), which conserve the total angular 
momentum in the symmetric gauge, AL Z = (An,- — Am,) = 0, will 
not be allowed. The only dynamical process allowed for this exciton 
(in first order and constrained to (i) and (ii)) is to decay in its own 
guiding center and create another exciton in a different one, as in (d). 
So, the number of particles in each guiding center m will be kept con- 
stant. The dispersion curves that we will obtain for the excitons are 
related to the breakdown of degeneracy between states like (a) and 
(d). Although the forward scattering approximation allows the de- 
cay processes like (b), we will not consider then due to an additional 
approximation made in the bosonization method. 



where N m = £~ = _ ro {c\ m c p .„ % - (Go| c\ m c^ m |G )) is the op- 
erator which counts the number of particles relative to the 
reference state \Gq) in the channel m. The last term gives 
the ground state energy (Gx|//o|Gn) of a system with a 
distribution of particles per channel given by the vector 
N = (No, ■■■,Nty-i) and the first one describes the excitations 
above this state. Formally, what allows us to rewrite (12.51 in 
this alternative bosonic representation, without a further ap- 
proximation, is the convenient linear dispersion relation in n 
of the system^. 

If we now include the Coulomb interaction term (12.61 ) in 
the forward scattering approximation and restrict ourselves to 
the spectrum of low energy excitations (close to the highest 
filled, V — 1, Landau level) the total Hamiltonian can be ex- 
pressed in terms of the bosonic operators bn) n , eq. (!3.H . Even 
restricting to processes that conserve the total number of elec- 
trons in each guiding center, as in the one dimensional case 25 , 
there is still more than one way to pair the fermionic oper- 
ators in (12.61 to yield the bosonic bilinear terms. If we set 
m\ = m\ and ni2 = m' 2 in (12.61 we get the direct scattering 
processes whereas exchange scattering processes are obtained 
for m\ = m'2 and m^ = m\ . This leads us to 



H i = -[ L [Sdi (n;m,m) + g ex (n;m,m)] b\ m b n ^ +H N , 

n>0.m.m' 

(3.4) 
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with 



g di [n; m,m) = n I dq$> m (q) (q) V£ (q) , (3.5) 

gex{n;m,m) = Ef d/ 8 m y (3.6) 
-n [ dq<t> m (q)<S> m i (q)V^(q) , 

-f/4 



where q = Iq, <J> m (q) = sfq~G m . m \\)e 



v n d M) 



E Self 



— z — |G v -i,v-i+n(q)| ■ 



dre ~?~L v _ l 



-V— 1+n 



(3.7) 



(3.8) 



dke-^L v - 1+n ( — )L w -i(—)jo(kq),(3-9) 



r = r/l and j7o (x) is the zero order Bessel function. 

In the above expressions, gdi(n;m,m') is the contribu- 
tion (see Appendix B for details) from the direct interaction 
process 4 ' 6 , that describes the recombination and destruction 
of the electron-hole pair, and g ex (n;m,m') is the contribu- 
tion from the exchange process (see Appendix B for details). 
Moreover, g ex (n;m,m') describes the bound state of a hole in 
the level v — 1 and an electron in the level V — 1 + n moving 
along different guiding centers 4 . This is the quantum ana- 
logue of the classical interaction of two particles of opposite 
charge under a magnetic field, which move parallel to one an- 
other with a constant velocity perpendicular to their separa- 
tion. This term has a constant diagonal contribution Eg el ^ that 
is the difference between the exchange self-energy of an elec- 
tron in the in the excited level v — 1 + n and the level v — 1 
from which it was removed. Besides, we see that the forward 
scattering approximation does not generate a diagonal inter- 
action in the guiding centers m. 

The term Hn in (13.4b contains only terms dependent on the 
number operator. Since they commute with the bosonic op- 
erators they have no influence on the system's dynamics but 
only on its chemical potential (see section [VII Al l. 



C. Diagonalization of the total Hamiltonian 

The total Hamiltonian in the forward scattering approxima- 
tion, 



Hj — Hq + Hi 



(3.10) 



«>0 m 



+ H N ,, 

n>0m,m' 



with g (n;m,m') = gdi (n;m,m r ) +g ex (n;m,m'), can be diago- 
nalized by an appropriate linear combination of guiding center 



bosonic operators b\ m . It is easy to see that such a combina- 
tion is given by 



<&m (q) 



qLm ( y 



,-? 2 /4 



bl(q)= £ *m(q)bl 
m=0 

(3.11) 

This new representation describes the degeneracy breaking 
between excitations like (a) and (d) in Fig.©. By using the 
completeness and orthogonality relations of the Q m {q) func- 
tions, 



m=0 



JO 



dq<£> m (q) (q) 



5{q 



(3.12) 



we see that the operators b n (q) also obey bosonic commuta- 
tions relations 



t)]=o, 

(3.13) 
(3-14) 



The same holds for the operator N (q), which plays the role of 
a "zero-mode boson": 



' y n (q) , b n , (q) J = §„,„'§ (q - q') , [b n (q) . 
and are related to the b\ m bosons by 

bt m = [ dq<t> m (q)bl(q). 



N(q)= £ <S>m(q)N n 

m=0 



(3.15) 



The total Hamiltonian can then be directly rewritten in this 
new representation as 



Hi 



n>0 



dqE n (q) b\ (q) b n (q) + H } 



N i 



(3.16) 



(3.17) 



with the energy spectrum 

E n (q) = n [h® c + V e " x (q)+n (q)) + E n Self . 

The fact that the Hamiltonian is diagonal in the basis b n (q) 
is a consequence of the forward scattering approximation. 
In the symmetric gauge, the system's angular momentum is 
given by L = (n, — mi) and is a constant of motion. There- 
fore, in (12.6b . the processes with non-zero matrix elements are 
the ones at which 

(n\ -mi) + (n 2 -m 2 ) = (n\ - mi) + (n 2 - m 2 ) . (3.18) 

In the forward scattering approximation we have m\ =m\ and 
m 2 = tn'2 ( or m i = m 2 anc ^ m 2 = rn[, in the exchange pro- 
cesses). Eq. d3.18b is then reduced to An\ + An 2 = 0, with 
Arii = n'i — rij. So, we are restricted to the first order terms 
in a degenerate perturbation theory, which present only terms 
that conserve the energy associated with the free portion, Hq, 
of the total Hamiltonian. In the bosonic representation this 
means terms like b\ (q)b n (q), but neither like b\ {q)b n i (q), 

with n 7^ n', nor like bn (q)t ' V (q). Formally this form ap- 
pears in our expressions after the angular integration over 
arg (q) in the matrix elements (see Appendix B), which clearly 
imposes the angular momentum conservation. 
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12 3 4 
\AL~=An-Am = o\ 



Figure 2: (a) Illustration of the linear momentum conservation by the 
action of b\ + b\ _ in the Tomonaga-Luttinger model, (b) Illustration 

of the angular momentum conservation by the action of the b\ m b n m * 
in our model. The particles are kept in their channels, therefore Am — 
0, and the variation on n , An, must be equal to zero too. 



IV. 



COMPARISON WITH THE TOMONAGA-LUTTINGER 
MODEL AND THE 2D-BOSONIZATION 



We can compare our Hamiltonian ( 13.101 ) with the one of the 
Tomonaga-Luttinger mode l 27 ' 28 : 



n>0 [j=± 

+ £ 82 (kn) (bl jbl f + b n jb nJ ^ 1 +H N , (4.1) 

associated with one-dimensional (ID) fermion systems. In 
this model the system's neutral excitations are distributed 
in two momentum channels, j = ± (left and right-moving 
fermions), presenting energies given by the linearization of 
the dispersion relation close to the two Fermi points, ±kf (see 
Fig-©)- In our two dimensional model the neutral excitations 
are distributed among the m = N§ guiding centers. Comparing 
(13.101 and (14.11 ) we see that, despite the different number of 
channels, we have the equivalence relations: 



g 4 (k n ) -» g(n;m,m) 



:(k n ) -^g{n;i 



In the Tomonaga-Luttinger model the terms associated with 
the interaction of particles in different channels, described by 
gi {k n ), present the form b\ -b^ n ■, or b n jb n ji. This is because 
of the linear momentum conservation in the system: b n+ cre- 
ates an excitation with wave vector kf + k n , while b\ _ creates 
an excitation with wave vector —kf — k n . So, b\ ,b ■, creates 
an excitation with total momentum equal to zero (see Fig.©). 

In the Tomonaga-Luttinger model this portion of the Hamil- 
tonian, associated with gi(k n ), does not commute with the 
kinetic energy term Hq = hvF'£n>oT.j=±k n bl jb n j. There- 
fore, it modifies the system's ground state, which has a dis- 
tribution of electron-hole pairs. The terms associated with 
g4 (k n ) describe the interaction between particles in the same 
channel. These terms commute with Hq and their influence 
is limited to removing degeneracies in the excitations^ 7 -. In 
our model, all the terms associated with the interparticle in- 
teraction, in the same or different channels, present the same 



form b* nm b nm i. As we have seen in the section above (after 
having made the forward scattering approximation) this form 
results from the angular momentum conservation, which does 
not allow the presence of terms like b\ m b , or b n , m b nm i (see 
Fig-©). Therefore, we will have a ground state like a Fermi 
sea, without particle-hole excitations or Landau levels mix- 
ing. All the interactions g (n;m,m'), with m = m! or m ^ m' , 
contribute to removing the degeneracies of the n levels. 

In a bosonized Hamiltonian, the presence of terms like b^b ' 
and bb, which separately do not conserve energy, "spreads" 
the single -particle spectral functions and gives rise to the 
anomalous dimension, a, in the correlations functions^. In 
the spin-full Tomonaga-Luttinger model, besides the spin- 
charge separation effect, the presence of this anomalous di- 
mension results in the absence of fermionic quasi-particles 
(vanishing of the quasi-particle residue Zk ~ \k — kf\ a as k — ► 
kf). In other words in the non-Fermi liquid behavior of the 
model. In the 2D bosonization of the electron gas 20 terms 
like b^b^ and bb, although present in the original interaction 
Hamiltonian, can be shown to vanish in the thermodynamic 
limit and in the asymptotic low-energy regime, giving rise to 
the Fermi liquid behavior. In our model it is clear that we will 
not have effects coming from an anomalous dimension, but 
we will have other effects that will give rise to a non-normal 
behavior. 



V. INCLUDING THE ELECTRON SPIN 

A. The Free Hamiltonian 

Including the electron spin, the free Hamiltonian (12.51 ) be- 
comes: 

#0 = faO C £ £ + V2)4,m/™- (5-1) 

s=± m =0 n=Q 

Considering the extension of the Hilbert space to negative en- 
ergy values, we can write Hq in terms of the new bosonic op- 
erators, b„ :S (q): 



Hq = H(O c £ / dq 



«>0 



K,s (q) = £ ®m (q) -t= £ 

m=l V n p=-o 



(5.2) 



(5.3) 



However, it is more convenient to define the following usual 
bosonic operators^ 7 -: 



1 
1 



bn{q) p = - 1 =[b n}+ (q) + b n -(q)} 



bn(q) a = -j=[ b n,+ {q)-b n -(q)] 



(5.4) 
(5.5) 



which are, as we will see, directly related to the charge and 
spin fluctuations in our system. Writing (15.2b in terms of these 
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operators, we have: 



separated: 



H = h(0 c £ / dq 



i=p,a 



Z j nbi(q) l b n (q) i +-N 2 (q) l 

n>0 L 



. (5.6) 



The bosonic operators b n (q) p and b„ (q) a commute. There- 
fore, they are associated with two different set of excitations. 
In one dimensional systems they describe the independent 
charge and spin fluctuation s 27 ' 29 ^ . In our system these op- 
erators play an equivalent role. If the higher Landau levels 
have both spin sublevels filled (see the next section) the sys- 
tem ground state will present total spin S = 0. Therefore, 
we will have triplet or singlet excitations. The triplet exci- 
tations, described by the operators b n (<?)p> gi ye f i se to poles 
in the density response function % p (q, to) whereas the singlet 

ones, described by the operators bn (<j) a , to poles in the spin- 
density response function % a (q, co) 4 -. 



B. The Zeeman Energy 

The inclusion of the electron spin brings together the 
Zeeman energy, H Ze = jVBgBZs=±Ln,mS c l,m,s c n,m,s, which 
splits each Landau level into two sublevels with energy sepa- 
ration of the order of AE Z = jp.BgB. However, we will treat 
only excitations between Landau levels. Our bosonic opera- 
tors b n (q) p or b n (q) a are associated with particle-hole pairs 
with the same spin, see eqs. (l5.3H5.5b . These are the only pos- 
sible excitations for a system with both spin sublevels filled 
(S = 0). The particle's spin will be conserved in the transi- 
tions and therefore the Zeeman energy will not have any in- 
fluence in the system dynamics. One can easily see that Hz e 
commutes with b n (q) p and b n {q) a . An appropriate descrip- 
tion of the intra-Landau levels excitations 31 , through an ap- 
propriate bosonization method, was also presented recently 
by one of us^ 2 -. 



C. The Interacting Hamiltonian 

We saw that the inclusion of the spin and the definition of 
the operators d5.4H5.5t bring the free Hamiltonian Ho ( 15.11 ) to 
the form (15.61 . Now, once the exchange processes occur only 
between particles with the same spin, the exchange Hamilto- 
nian Hf x , dB3b . will be diagonal in the s index, like Ho, (I5.ll ). 
Therefore, it will be transformed in the same way as Ho, be- 
ing written in terms of b n (q) p and b n {q) a . On the other hand, 
the direct interaction Hj, (IB 2b . will not be diagonal in the 
spin index. The sum over the spin indices, £ s y, results in 
the cancellation of the terms proportional to b n {q) a and in the 
presence of an interaction equal to 2VJ. (q) associated with 
the charge operators b n (q) . Physically, the direct interaction 
gives rise only to charge excitations, while the exchange in- 
teraction gives rise to charge and spin excitations. The final 
Hamiltonian will then present the charge and spin fluctuations 



H T = H {p) +H {a \ (5.7) 
H« = £ [dqE®(q)btiq) t b a (g)o 

n>Q J 



E n i] (q) = n[h® c + 2^ p V^ + V?M)+E n Se i f ,(5X) 

with i = p,a and the terms of the energy spectrum En (q) 
given in d3.7H3.9l ). We did not write the terms proportional to 
the number operators. In the limit ^->»we will have the 

asymptotic value En = nh(O c + Eg e ^ since V" x (q) and (q) 
goes to zero. 



VI. SPECTRUM OF NEUTRAL EXCITATIONS 

For the first excited level, n = 1 , the neutral excitation spec- 
trum that we have obtained , eq. d5.8l ), is identical to the one 
determined in eq.(4.2) of the Kallin and Halperin's paper 4 . We 
just have to perform the angular integration over arg (r) in the 
latter to obtain (15- 8b . For n > 1 the only difference is the n 
factor multiplying V7. (q) and V" x (q), that we will discuss be- 
low. The authors of Ref.[4] clarify that (15.8b . within the strong 
field approximation, applies for two cases, at which it is exact 
to first order in (e 2 /e/) /h(O c : (i) excitations associated with 
n = 1 in systems with arbitrary V and (ii) excitations associ- 
ated with arbitrary n since the system presents V = 1 ; in both 
cases the system must present the two spin Landau sublevels 
filled, as we suppose in sec. dVBb . Despite these restrictions, 
these are the excitations of major interest. Outside the do- 
main of these conditions eq. d5.8b represents a less controlled 
approximation, once the multi-exciton states were neglected. 

In the bosonization procedure we can easily understand the 
origin of these restrictions. To obtain a precise bosonic de- 
scription of the excitations we have approximated all the ma- 
trix elements G 



p,p-\ 



that describe the excitation of an elec- 
tron to the p + n level, leaving a hole in the p level, by the 
Fermi level matrix element, Gp,p+n\ p=v _i = Gv-i,v-i+n, see 
Appendix B. In this case we have a precise description for the 
excitations associated with n = 1 , as they are represented by 
the unique element G v -i.v> and for the arbitrary excitations in 
a system presenting V = 1 , since they are all associated with 
the same element Go. n - In these two cases all the particle-hole 
pairs that contribute to a given excitation are associated with 
the same Landau levels. So, we do not have to consider all the 
n Landau levels between v — n and V — 1 that would contribute 
to a given E„(q) if V — n > 4 . However, once we introduce 
negative energy states, even for filling V = 1 we will have a 
contribution of n states. This is the origin of the n factor mul- 
tiplying VI (q) and V" x (q) in (15. 8b - But, as we will see, the 
n = 1 magnetoplasma mode will be the most import excita- 
tion in our analysis and, in the strong field limit vfUO c 3> e 2 /el, 
the higher magnetoplasmon modes will not present a relevant 
modification due this n factor. 

We justify the approximation in the matrix elements, 
Gp, P +n\ p=v _i = G\-l,\-i+n, with the strong field approxima- 
tion, see Appendix B. However, we should mention that it is 
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also a good approximation to describe the low energy excita- 
tions (n«v) of a system in low magnetic field (v S> 1). In 
this case, for small wave vectors q, it can be shown that the 
matrix elements G v -i v-l+n (?) depend approximately only 
on n and the value adopted by V is irrelevant. This property 
was explored in a similar bosonization method 16 . 

The dispersion curves for the modes close to ftia c are ob- 
tained calculating ( 15.81 ) for n = 1 . For V = l 3 -, e.g., we got: 



We get a fermionic field appropriate for a bosonic represen- 
tation with b njn , introducing the phase representation 16 



Vm (9) = 



1 



(7.5) 



7(0 

J V=n=l 



(<?) 



8/2 V 2 



-tnt 



« 2 



/o 



25 iiP \l-qe t 



£ 

4 
(6.1) 



In this way the set of operators {c n . m ,b n , m ,N m ,\\f m (Q)} is 
completely analogous to the set {ck n j,b n j,Nj,\\ij(x)} if we 
make the replacements: j —>m,x—> 0, k„ —> n, ^-Nj —> N m 

Therefore, writing the fermionic field 



and 



/2k' 



\f m (0) in terms of the operators b n (q) of ( 13.1 11 ), it presents 
the following representation, analogous to (I7.2H7.4I ): 



These expressions, for the charge and spin modes, i = p,a, 
are exactly the eqs. (4.4) and (4.8) in Ref.[4] or (3.3) in 
Ref . [31]. The first term between curly brackets is independent 
of q and corresponds to E^j-. The second term corresponds 

to V^ =1 (q) and the last one to 25 (i pV^ (q). For higher n or 
V values we have much longer expressions, so we do not ex- 
plicitly show them here. 



VII. SINGLE-PARTICLE PROPERTIES 

A. Fermionic Operator 

The bosonic representation of fermionic fields, also known 
as the Mattis-Mandelstam formula, is an identity in one di- 
mensional systems if the theory can be formulated in terms of 

a set of fermionic operators, |o;j,cj[/^j = which 

are labeled by the "species" index j = l,....,Nj and a discrete 
and unbounded "momentum" index k, which typically labels 
the eingenenergies, E\, of the non-interacting system 18 . The 
fermionic fields, (x), are defined as: 



ft 



e~ iKx c 



n,Ji 



(7.1) 



where L is the length associated with the system size and c„j 
remove the fermion from the state k n = ^n. They have the 
following bosonic representation 18 ^: 



1 



-ek„/2 



,-ik„x } 



k„>0 V 
1 



n,j I ' 



,t 



(7.2) 
(7.3) 
(7.4) 



where b n j are bosonic operators 18 and £ > is an infinitesimal 
regularization parameter such that £ — > + at the end of the 
calculations. 



M>m (9) 
©m (9) 



1 



lim ■ 



-UmC e , 



i / dq<t>,„ (q) Y 

J o „=1 



-me 



b n (q) 



(7.6) 



(7.7) 



In order to determine the time evolution of the fermion field 
( 1X61 ) we first note that b n (q,t) = e~ mn ^'b n (q), with E„ (q) 
given by d3.17| l, and that the time evolution of U m is given by 



U m (t) (/,„<< " 



' m t i(nZi-H%\v ,,W, 

N l — e \ N N W m +lJ JJ 



(7.8) 



where H™ is the portion of the total Hamiltonian propor- 
tional to the N m operator. Acting on the ground state we get 
(0| U m (t) = e^' 1 - 1 " 1 ' (0| U m , where ^ corresponds to the chem- 
ical potential associated with the creation of a particle in the 
mth channel, ,u+ = Eq (N m + 1) — Eq (N m ). Grouping the terms 
proportional to N in (13.101 . which are given in the Appendix 
B, we get 

rt = toc -I VM )) + ^^W. (7.9) 



B. Retarded Green's Functions at T — 

The retarded Green's function, at temperature T = 0, as- 
sociated with a particle in the quantum state \n,m) (with n 
measured from the Fermi level, see Fig.Q) is given by^: 



iG n% ( f ) 



*Ve(0 
iG<Jt) 



(0)+4,m (0) 



/•2n 

10(0 / dOe 
■JO 

•t 



m9 



G>At) + G< (t) 



,(7.10) 



v«(9,Oy2,(o,o) 

(0,0)^ (9,0 



where we use the relation (17.51 between the fields \|/ m (0) 
and c„, m . The calculation of G< Q (t), from the representation 
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( 17.61 ), is presented in Appendix C. We get: 



271 



iGfLit) = 9(0 / 

JO 



(7.11) 



* me( = L e ^r^- M i^^)e-^\ (712) 

271 

where corresponds to the chemical potential associated 
with the removal of a system's particle. 

We can calculate the retarded Green's function associated 
with different states. For a particle in a position r we have 



iG Ret (r,t) = 0(O(v(r,O¥ t (r,O) + ^ t (r,O) ¥ (r,f 

1 2 r'Ret 



n'=0 m=0 



(7.13) 



where we use the representation of \|/ (r) in the Landau level 
basis, ( 12.21 ), and (r|«',/n) is given by (I2.3I ). Notice that we 
defined (f), eq. (17.1 II ). with the index n relative to the 
Fermi level, according to Figs.© and ©, whereas the index 
n' above is relative to the single particle ground state. So we 
must write n' = n — (v — 1) in G R f m (t), getting 

G Re '(r,t) = tt \(r\n,m}\ 2 G R n ^_ 1); Jt). (7.14) 

n=0 m=0 

If we want the Green's function associated with a particle 
in a different state, e.g., the wave vector q, we have to replace 
(r\n,m) by (q\n,m). 



C. Single-Particle Spectral Functions 

Now we will determine the single-particle spectral func- 
tions, A (go), which gives the spectral distribution generated 
by the inclusion or removal of a particle in the system. The 
spectral function associated with a particle in the sate \n,m) is 
given by 

1 1 r°° 

A n . m (©) = - -lmG% (CO) = —Im dte~ im G R % (f ) . 

(7.15) 

We need the Fourier transform of the function G Ret m (t). We 
then take n > 0, creating a particle or a hole above the Fermi 
level (see Fig.©). Since our ground state does not present 
electrons above the Fermi level, see section dlVI ), only the 
first term between square brackets of (17.111 ) will contribute to 
A n >o,m (to), describing the creation of an electron in the sys- 
tem. Substituting the expressions (17.1 II ) and (17.121 in (17.151 ). 
we get 

1 f 2lZ 4- 



Now, the exponential function of the sum X^° =1 can be ex- 
panded in a series that, after the angular integration over 9, 
will result in a finite number of terms for each n. It is easy to 
see that 



l(J n=l,m 


(0 


= B(t)e 


-ip + t 

i 






(7.16) 


■ r Ret 


(0 


= B(t)e 


-i M +t r 

Jo 


dq<s>l(q)e 


iEi(q)t 

? 


(7.17) 


iG Re \ 


(0 


= \m 


e ~ ¥+t J 


[ dq<Pl(q) 



~ e -iE 2 {q)t 








dq<S>l 


(q') e -'[ E ^h 







and so on for n > 3. Therefore, we see that for each iG Re j n (t) 
we only have the contribution of the normal modes presenting 
n < n. This behavior is general, occurring in any bosonized 
system that presents a discrete levels structure. Recently, it 
has also been observed in the treatment of the Luttinger liquid 
in a finite one-dimensional wire 36 , where the discrete levels 
structure came from the box-like boundary conditions. We 
can also see that we do not have the contribution of the mode 
h = n. Therefore, according to (17.151 ) and (17.161 . the spectral 
function associated with the creation of a particle in the first 
excited state, n = 1, is given simply by 



A„=i, m (co) = 8(co-ai h 



(7.1 J 



In this situation the only energy level excited by the intro- 
duced particle is the one given by the chemical potential, with- 
out any contribution from the excitation of bosonic modes. 
This was expected as a consequence of the first order ap- 
proximation. Consider the inclusion of a particle in a state 
\n = l,m) , as presented in Fig.((3]a). Again, since our ground 
state has no particle-hole pairs, there is no empty state to 
which this particle could relax. Therefore, it cannot excite (or 
decay into) any of the bosonic modes of the system. The par- 
ticle can neither move between the others channels, once the 
number of particles in each channel must be conserved. These 
means that the particle can not "see" the degeneracy breaking 
in m, associated with the dispersion AE n= \ (q) in the level it is 
included. It is easy to see that this is as a consequence of the 
first order perturbation theory. Due to the angular momentum 
conservation, the particle could jump to a channel m' ^ m, e.g. 
m' = m + 1 , only exciting another particle in another level n 
above it, as presented in Fig.©. However, this should be a 
virtual transition, once it does not conserve energy. 

According (17.151 ) and (17.17b . the other spectral functions 
will be given by: 



J £/$*£($) 5 (o> -El® -/!+), 

00 dq® 2 m (q) { [5 (CO - E 2 (q) ^+) + (7.19) 
dq'& m {if) 8 (CO - Ex (q) E X (if) » + ) 



1 

2 Jo 



xe 



A n =2,m (CO) 
A„=3,m (CO) 



and so on. It is clear, from the scheme of Fig.®, that the 
above expressions correspond to the distribution (or decay) of 
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Figure 3: Scheme representing the distribution, in the bosonic 
modes, of the energy of the particle introduced in the system: (a) 
a particle created in the n = 1 level can not excite any bosonic mode, 
(b) created in the level n = 2 can excite the bosonic modes associated 
with n = 1, E n= i (q), (c) created in the n — 3 can excite the bosonic 
modes associated with n = 1, £'„ = i (<?), and n = 2, £,,=2 (<?), and so 
on. 




2 An = 1 




functions (see Fig.©). 

Analyzing the above expansion of G Ret m (?) and recalling 
that Jq dq<&} n (q) = 1 we see that A njn (co) obeys the sum rule 
f^dtoA^m (co) = 1, for all {n,m}. 

The spectral density associated with the inclusion of a par- 
ticle on the position r is given by 

A(r,(0) = --ImG Re ' (r,a) = £ £ |(r|?i,m)| 2 A„_ (v _ 1)!ffl (co 



n=0 m=0 



(7.20) 

where we use the relation ( 17.141 ) between G Ret (r,t) and 

G Mv-i),m(0- We verify that /-»^»A(r,C0) = 1, once we 
have the correct normalization £„ m |(r|?i,m)| 2 = 1. Particu- 
larly, for the creation of a particle at the origin, r = 0, we have 



A (r = 0,0)) 



1 



(271/) „ =0 



Y, A n-(v-l),n 



CO 



(7.21) 



because (0\n,m) = 8„. m / V2nl 2 . One particle created at r = 
add to the system an angular momentum AL Z = 0. There- 
fore, it must decay only into the modes with n = m, which 
also present zero angular momentum and are indeed the only 
terms left in (17.20b . We will see the graphic representation of 
A (r = 0,co) in Fig.©, after the inclusion of the spin degree of 
freedom. 



D. Including the Electron Spin 



Figure 4: Scheme representing the "difficulty" (impossibility in first 
order) of a particle added to the system in the first excited level (a) 
to make a transition to the neighbor channel, with Am — +1 (b). To 
conserve angular momentum, AL Z = An — Am = 0, we must have 
another particle making a transition An = +1 simultaneously. This 
must be a virtual transition, once it does not conserve the energy. 



a single-particle excitation in a set of collective bosonic exci- 
tations of total energy close to the free energy of the particle 
introduced. As the energy of the introduced particle increases 
it excites a larger number of collective modes. This kind of 
distribution is the essence of the bosonization methods^. In 
the next section, after we include the particle spin, we will 
present the graphic representation of some of this spectral 



The details of the introduction of the electron spin, s = ±1, 
are presented in the Appendix D. For the Green's function we 
get 



iG n%,s ( f ) 



M„ 



r2% 

9(0 / dQe 
Jo 



/«9 



-e-^'MU(t) 



(7.22) 



1 2l,T=l ^- i " 9 /0°°^,(9)^-'' £ " P)( * ) '+«-'' £ " O)( * ) ') 

2n e 



with En given by d5.8t . 

The spectral function A nmiS (co) presents the same structure 
of the equations (l7.18H7.19l under the following substitutions: 



8((0-E n (q)-^ 
5((0-E n (q)-E n (q')-^ 



i=p,a 



(0 



- £ 8(co-4 iJ (?)-M 



(7.23) 



£ 5 U-E n '> (q)-E [ '> (q') -4) +25 (w-E^ (q)-E^ (q') 
/=p,a 



and so on. The functions A nms (co) also obey the sum rule J^ddtAn^^ (co) = 1, for all {n,m,s}. 

In Fig.© we present examples of the form acquired by these spectral functions. The spectral function A„ = i , m i (co) corresponds 
simply to a 5— function centered at the chemical potential energy /uf. From A„ = 2, m ,. ? (co) on to higher values of n we start to have a 
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complex energy distribution. The particle excites simultaneously, according the scheme of Fig.([3]b), the charge and spin modes, 
associated with (q) and E^2\ The spectral function A n =2,m,s ( w ) ( see (17. 19b ) can also be written as: 



where we have performed the integration J °° dq = J* °° 



and 



| ^(9) | 

means a infinite sum over all the q val- 



^n=l,m=0,s( m> 



v=2 B=5T 



-fr-^™! (§)+«+ = 

ues such that (<?) +^ corresponds to the given 0). So 
A n =2,m,j wm present peaks (integrable divergences), with 

weights proportional to <J>^ (q) /2 where dE^l^ (q) /dq = 
and the system presents maxima in its density of states. This 
behavior is similar to the behavior of the dynamical structure 
factor, 5 (q, 0)), measured by inelastic light-scattering 710 in in- 
teger quantum Hall system, which confirm the existence of 
the rotons' minima in the dispersion of the collective modes. 
We can easily calculate S(q, CO) within the bosonization pro- 
cedure, but since we obtain the same results presented in 7 we 
do not present the calculation here. 

The expression (17.241 ) is useful to understand the behav- 
ior of A njn s (to). However, the more direct way to evaluate 

it is substituting 5 ^00 — ^n=i (*?) — Ms" J by a Lorentzian of 
negligible width. This procedure allows us to compare the 
weight of the different peaks, while keeping the area under 
each one of them constant. In Fig.© we present the behavior 
of A n= 2,m=o,s (o) for two different filling factors V. We see 
peaks in the energies of the extreme points of the charge and 
spin fluctuation modes and the energy of the particle being 

distributed between the minimum of (q) and the maxi- 
mum of (<?). This behavior is analogous to the distribu- 
tion of energy of an added particle in the Tomonaga-Luttinger 
model 37 . There the energy of a particle with momentum q is 
partially distributed between the charge and spin modes, pre- 
senting dispersions v p q and v a q, and also over energies above 
and below these ones, due to the existence of the anomalous 
dimension. Here we do not have any effect of an anomalous 
dimension, see sec. dlVt . however, both dispersion relations, 

E^ (q) and E^ (q), present a non-monotonic behavior, giv- 
ing rise to the many peaks in the spectral function. We also 
have the presence of the function <E>^, (q), which distributes the 
energy according to the relation between an excitation in the 
guiding center m and one with wave vector q, which we will 
discuss in see l VII Fl 



A n=2,m=0,s(<°) 



(w-i4)/a> c 



0.5 1 1.5 2 2.5 



(W-M*)/(0 C 



(m-f4)/m c 

Figure 5: Examples of spectral functions A nm= Q S (co) for a system 
with filling factor v = 2 and a typical transverse magnetic field of 
5T in a GaAs heterojunction, where we have h(£> c [K] ~ 20B [T] and 
e 2 /el [K] ~ 50 y/B [T]. The 8-functions in flT23) where substituted 
by a Lorentzians of width 10 _3 O) c . 



where /Uq = hx£> c + sAE z , AE Z is the Zeeman energy, and the 
retarded Green's function would be given by 



i 



co-(»-l)w c -^ + /ri 



(7.25) 



Now let us study what happens in the interacting system. For 
n = 1, the introduction of the interaction results in a simple 
renormalization of the chemical potential /Vq to fif, similar 
to d7.9t . For n > 1 we have the spectral functions with the 
(integrable) divergences of Figs.© and (J6J). If we take the 
Fourier transform of G£ff 2 m . it), given by (17.221 ). we get 



G n=2,m,s ( W ) 



2 u 

t=P,o 



dq *m . 

" a>-E^ =l (q)-^f + iy] 



E. Non-normal System 

In a non-interacting system, for n > 1, we would have 
(»=1) 

(co) = 5 (co - (n - 1 ) co c - , 



So, instead of a single pole, as in (17.251 . we have a continuum 
of poles which give rise to some integrable divergences in the 
spectral functions. Therefore, we see that the effect of the in- 
teraction on the non-interacting spectral function is very dif- 
ferent from shifting and broadening the non-interacting peak. 
It generates peaks with different energies and weights, as- 
sociated with both charge and spin fluctuations. In normal 
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{<o-a> e -i4)/e 2 M 



(<a-a> c -Ms)/eV£l 



Figure 6: Spectral functions A n= 2 m =o s (<*>) f° r different filling fac- 
tors v. In the second line the 8-functions in ( I7.23t where substituted 
by Lorentzians of width 10 _6 CO C , stressing the integrable divergence 

at the local extremes of the dispersions (q). The charge (contin- 
uous) and spin (dotted) contributions are represented separately. The 
third line represents the addition of the two previous contributions, 
but with Lorentzians of width 10 _3 C0 C , such that we can see the rela- 
tive weight of the different peaks. Horizontal and vertical lines with 
the same dotted pattern correspond to the same energy (co = E with 
h=l). 



systems 25 ' 38 , for each non-interacting eigenstate, \n,m,s), we 
would obtain a single solution of the equation 



CO - (n : - 1 ) C0 C - fjf - ReL re! (\n,m,s) ,C0 - fj' s 



0. 



and this would be the energy of the quasi-particle in the inter- 
acting system. Here E ref is the retarded self-energy 25 ' 38 . How- 
ever, in our case the above equation presents many solutions, 
one for each divergent peak in A„. m s (co). So, in general, we 
see that we have a non-normal system. However, in our treat- 
ment, there are two exceptions. The first is the case associated 
with the inclusion of a particle in the lower allowed energy, 
when A n= i jm ^(co) is given by a simple 5-function. The sec- 
ond is the case A n !m -,N<s>->l,s that we will see in the next 
section. 



F. Influence of the Guiding Center Position 

So far, in analyzing the spectral functions we have concen- 
trated on the influence of the varying densities of states. How- 



ever, we have also mentioned the influence of the guiding cen- 
ter position, defined by m. The function <E> 2 (q) weights the 
importance of the density of states at each point. The spectral 
functions presented in Figs.© and ([6} were calculated with 
m = 0. In this case <I> 2 _ (q) does not present oscillations and 
limits the influence of the densities of states to the range of 
medium (not too large) q s. However, once we increase the 
value of m, the function <J> 2 (q) spreads and starts to present 
m + 1 oscillations, with the last maximum being the largest 
one, see Fig.©. Once the dispersion relations became softer 
at large qs this region starts to concentrate most part of the 
spectral weight. Consequently, as shown in Fig.(JTJ, the spec- 
tral function A„ = 2. m ..s (co) starts to present weight at larger fre- 
quencies, with a peak much more pronounced that the one in 
<E> 2 (q). In the limit m — > N® S> 1 the peak will be concen- 
trated at the asymptotic value C0 C +E^J, with a char- 
acteristic modulation. It is clear that, in this limit, more gen- 
eral spectral functions, A„. m s (co), would also present them- 
selves concentrated at well defined energies (shifted from the 
non-interacting values) and we will not see any effect of the 
decomposition into charge and spin fluctuations, once both 
present the same asymptotic limit. Systems presenting higher 
filling factors, as in Fig.©, present a larger number of peaks, 

because of the V oscillations in (q) and (q), but their 
dispersions relations still converge to their asymptotic values 
in the limit m — > N® 2> 1 - 

We believe that this result can be understood as follows. A 
particle in a state \n,m} presents a probability density of width 
\/2n + 1/, spread over a circle of radius \/2m + l/ 39 . There- 
fore, for not too large values of n, creating a particle in a state 
written in terms of m~ N<t> means to create a particle in a large 
radius, close to the edge of the system. Most of the electron- 
hole pairs excited by the introduction of this particle should 
also spread over large radia, where their net wave function 
overlap is small and they do not fell much the presence of 
each other. Therefore, the energies excited in the system will 
be associated with electrons and holes far apart, E$ e if- How- 
ever, the particles created close to the system's origin, with 
states written in terms of small m's, would present a much 
more complex structure in their spectral functions, because the 
electron-hole pairs excited will present a much larger correla- 
tion. Below, we will explicitly calculate the spectral function 
of a particle added at the system's origin. 

The inclusion of the spin in G Re< (r,t) is made by substitut- 



ing Gg£(0 by in 



--Ret 



*) = -»Tl 



,r \n,m 



\ 2 G Re ' 



(v— l),m,s 



The same must be done for the calculation of A (r,s,co). For 
a particle added at the system's origin we have the expression 
equivalent to (I7.2U : 



A (r = 0,s, co) = — £ A„_ (v _! w (co) . 

(271/) n=0 



(7.26) 



Therefore, A(r = 0,s,co) will present many peaks. This is 
similar to the sum of the different peaks presented in Fig.©, 
but with each peak being modulated by the corresponding 
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A n=2,m,s<a>> 



A n=2,m,S<U» 




0.2 D.4 0.6 



i A. 




q (lo-fOr-ftl/eVel 

Figure 7: Spectral function A n —% m (<*>) f° r m — 0,3 and 15, in a sys- 
tem presenting a filling factor V = 1. Horizontal and vertical lines 
with the same dotted pattern correspond to the same energy. The 
vertical line indicated as Eg e i f would be the position of the peak in 
the limit m — > N& > 1 . 



<E>m = „ (q) instead of the same O^ =0 (q). The terms presenting 
n < V — 1 will generate contributions to frequencies 0) < 0, 
describing the effect of creating a hole in the system. They 
are similar to the ones in the region CO > 0. Fig.© presents 
the first two peaks of A (r = 0,5, Co) in the region CO > as an 
example. 

Expression d7.26t also indicates that as the filling factor in- 
creases the spectral function develops better defined peaks in 
the CO > region. In a system with V 3> 1 the spectral function 
A (r = 0,5, CO) at frequencies CO > C0 C (associated with \n\ > V) 
will be governed by the functions ^„=n>v>i (<?)■ So the sys- 
tem will tend to present a set of well defined peaks, shifted 
from the non-interacting values, as in a normal system. 



VIII. CONCLUSIONS 
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Figure 8: Spectral functions A n= 2ms (co) for different m values in 
a system presenting filling factors v — 2 and v = 6. The vertical 
dotted lines correspond to local extremes of the normal modes, as 
in the figure before, except for the vertical line indicated as E$eif, 
corresponding to the asymptotic value. 



(2nl) A(r=0,s,O)) 




Figure 9: Spectral function for a particle created at the center of a 
GaAs he teroj unction presenting filling factor v = 2 under a trans- 
verse magnetic field B — IT. 



In this work we presented a model and a method for the 
study of the integer quantum Hall systems presenting both 
spin sublevels filled (S = 0). We showed that, to first order in 
y = (e 2 /e/) / h(ti c , our model of interacting one-dimensional 
channels within the forward scattering approximation, which 
keeps constant the net number of particles in each channel, 
is equivalent to the strong field approximation. We saw how 
this model can be solved, with some well controlled approxi- 
mations, through an appropriate variation of the Landau level 
bosonization method. Using this set of model and method we 
determined the collective excitations of the system: spin and 



charge fluctuations, dependent on the direct and exchange in- 
teractions. We obtained the same dispersion relations as those 
presented by Kallin and Halperin 4 and MacDonald 5 for the 
excitations close to the cyclotron frequency C0 C . 

From these non-monotonic dispersion relations and vari- 
able density of states, we determine the single particle spectral 
functions through the bosonic Mattis-Mandelstam represen- 
tation of the fermionic operator. For states created close to 
the central region of the system, these functions presented the 
characteristic behavior of a strongly correlated non-normal 
system, including a spin and charge spectral decomposition 
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and absence of quasi-particle excitations that keep similarities 
with the ones of a non-interacting system. However, for states 
close to the edge of the system these functions tend to present 
single peaks, shifted from the non-interacting values and with 
a characteristic envelope, like in a normal system. It is clear 
that, in our treatment, we do not consider the effects of the 
confining potential in the Hall system, which generate gapless 
edge excitations 40 and are fundamental to a complete descrip- 
tion of the edge states . However, independently of the addi- 
tional effects of the confining potential, what we see is that the 
bulk quasi-particle states seem to be much more strongly cou- 
pled to the collective excitations in the system than the edge 
ones, which seems to behave like normal Fermi liquid parti- 
cles, even in higher energies. 

The experimental verification of the results we predicted 
seems to require a control over the state of the created par- 
ticles higher than the one available in current experiments. 
However, if as we predicted here, the spectral functions for 
states created at the central region and those close to edge of 
the system have such a distinct behavior, this should not be 
too difficult to observe. 

As we mentioned in the introduction and in Sec.(III.A), 
one of the advantages of the present model for the quan- 
tum Hall systems is that it presents a prescription for its im- 
provement, once there are methods developed to treat the 
net particle transfer between neighbor channels in a bosonic 
description 41 . We saw that the previous perturbative methods 
for IQH systems 4 - 5 used y as the expansion parameter. How- 
ever, they do not give a systematic description, order by order, 
in this parameter. Indeed, not even the first order in y is fully 
described, indicating that this may not be the better expansion 
parameter. So, we propose, start from the forward scattering 
approximation, which is already a good initial description as 
we saw, and try to get higher order corrections taking the net 
particle transfer between "m" neighbor channels as the expan- 
sion parameter. Once this transfer would be between channels 
with non-zero relative angular momentum, we suspect that the 
effective transfer potential can be related with some kind of 
Haldane pseudo-potential 42 , which, at least in the fractional 
quantum Hall systems, are known to be the really effective ex- 
pansion parameters. Certainly, with the inclusion of the trans- 
fer between neighbor channels, we will have a Hamiltonian 
whose ground state will present particle-holes pairs, taking 
into account the Landau level mixing, which would be very 
convenient in many applications. In this case, the spectral 
functions A n= \ (co) should already present a complex struc- 
ture, in contrast to the 5-function that we have got. 

We also mentioned the fundamental importance of the col- 
lective modes in the FQHE, evidenced by the Hamiltonian 
theory of Shankar and Murthy 11 , which have many similar- 
ities with a bosonization method. Their fermionic field as- 
sociated with the composite particle, \|/cp> is very similar to a 
Mattis-Mandelstam representation associated with the magne- 
toplasmon excitation. This seems to us another evidence that 
systematic bosonization methods created for quantum Hall 
system should be explored. 

Finally, we note that the one dimensional characteristics 
that seem to be retained in the 2DEG under a strong perpen- 



dicular magnetic field, the influence of the collective modes in 
the quantum Hall quasi-particle properties and the possibility 
of using the relative angular momentum between particles as 
an expansion parameter in a field theory, probably can still be 
further explored in quantum Hall systems. Here we presented 
a method that partially explores the first two characteristics, 
indicates a way that can be followed in the investigation of the 
third one and can be improved for further applications. 
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Appendix A: SOME PROPERTIES OF THE FUNCTIONS 

G n , m ' (q) 

In this Appendix we just list some relations among the 
Gn,m (q) functions 21 , with q = q x + iq y , which are useful for 
the matrix element manipulation. 

Hermitian Conjugate: 

G m ', m (q) = G* m m , (-q) = G m . m , (q*) . (Al) 
Matrix Product: 

^G m / i/ (q)G /)OT (q') = e^G^q + q'). (A2) 
Full Landau Level Rule: 

£ G mm (q) = A/*8 q , = ^-#*8 (q) . (A3) 

m=0 * 

Orthogonality: 

J d 2 qe~ q2/2 G mjl1 > (q)G„y (q*) = — 8 m , n 8 m /y. 
Matrix Elements: 



(n,m\e ,qr \n' ,m') = e ^G„y (Zq*)G m y (/q) . (A4) 

Appendix B: THE INTERACTION MATRIX ELEMENTS 

To obtain the matrix elements of the interaction Hamilto- 
nian in the forward scattering approximation, eq. (13.4b . one 
has to pair the creation and annihilation fermionic operators 
of (12.6b in all possible ways and set the same guiding center 
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m for both. This is easily accomplished by defining a general- 
ized density operator 



p(q,x) = Hr £ G„y(Q*)G„ 

{«,m} 



(q)4,m c »'m', (Bl) 



with q = q/ , Q = q + z'x and x = x/7, and finding its bosonic 
representation in terms of the operators b n>m or even directly 
in terms of b n (q), as we will do here. It is straightforward, 
though a little cumbersome, to show that the two ways of pair- 
ing the creation and annihilation operators of the Hamiltonian 
that yield the direct and exchange terms as follows 



H, 



Dir 



1 1 



2 (2nY 



d 2 qV(q)pHq)p(q)-±-V(0)N T , (B2) 



H 



Ex 



1 1 



J 2 



d 2 q / d 2 xV (i)«-T V'(q-**-q*-S)/2 



1 



xp(q,x)p T (q,x) + -y(0)^, 



(B3) 



where p (q) = p(q,x = 0) is the electron density operator 
and we introduced the total particle number operator Nt = 



Hn,m c n,m c n,m- To obtain the exchange term , the matrix ele- 
ment in ( 12.61 ) must be expressed in the position basis, |r,r'), 
using (12.31 . Then, in the resulting expression, we replace the 
integration variables, r and r', by the relative, x = r — r', and 
center of mass, X = (r + r') /2, coordinates evaluate the inte- 
gral with respect to the latter and use the relation (lA2b from 
Appendix A. 

To obtain the bosonic representation of p(q,x) we first 
rewrite the sum over {n,n'} as £„„/ = £„>„/ +L«<«' +(« = 



n'), doing n = p + n and n' = n in the terms where n > n' 
and n~^n' when n < n'. At the same time we can make the 
forward scattering approximation, m = m', to get 



p(q,x) = e t £ G P)P+ „(Q*)£G OT;W (q) cJ m Cp +niOT 
n>0,p L m 

+G P+ n,p (Q*) £ G m , m (q) 4 + „ !m c p , m (B4) 

m 

+Pw (q,x), 

where 

p w (q,x) = e'^ Y^G^n (Q*) G mm (q) 4,m c n,m • 



Now we approximate G PiP+n (Q*) ~ G P)P+ „ (Q*) | p=v _ 1 
since we are restricted to the low energy excitations close to 
the higher filled Landau level, V — 1 . Using simultaneously 
(|3~TT ) and (|3~TT1 i. we get: 



p(q,x) 



« 2 



^£v^[Gv-l,v-l+» (B5) 
v<7 « 



+ Pn (q,x). 



+G v _ 1+n , v _ 1 (Q*)^(?) 
The bosonization of pyv (q, x) is done in the same way as the 
bosonization of Hqt^, eq. (13.31 l. because both have the same 
form except for the expression G„,„ (Q) G mjn (q) which re- 
places the simpler function n. Calculating the commutator 



[Pn (q,x) ,b NM ] = Gmm (q) —?=■ Y I g p,p (Q*) ~ G p +n, p +n (Q*)] 4,m c p 



we finally conclude that pyv (q,x) must have the following bosonic representation: 



p N (q,x) = (G 1 p N (q,x) |G ) +e tt £ ( i- 



u n.m u n,nn 



(B6) 



where 



v-l ^ / r 2 \ 

(G | p w (q,x) |G ) = £ G„,„ (Zq* - ix*/l) £ G m , m (Zq) = -^L^ -^ 5 

n=0 m=0 ' \ Zi / 



(q) 



In the last step we have used the relations (IA3t and J^_o^ (•*) = ( x ) • 

The bosonic representation of p (q) is a particular case of relation (IB lb and is given by 



p(q) 



_ Y V" G v -i. v -i +n (q*)b n (q) + G v -i +ny -i(q*)bl(q) + p N (q). 
v<7 „>o L J 



(B7) 



(B8) 



15 



The last term, pyy (q) = Pn (q^O)* wm gi ve a nui l contribution after its substitution in (IB 2b if we neglect cubic terms in the boson 
operators. These would describe the dissociation of an electron-hole pair—. 

Substituting dB8b in (IB 2b . neglecting the cubic terms, considering the system neutrality (that cancels the effect of the mean 
electron density (Go| Pn (q) |Go) = ^5(q) present in p^ (q) = p^ (q>0), eq. (lB6b ) and performing the angular integration over 
arg (q) (that eliminates the non-diagonal quadratic terms, see section UlI Cb . we finally obtain the direct term of (13.41 . 

Now, following the same procedure and substituting the bosonic representation, dB5l > and dB6l ), of p (q,x) in (lB3b we verify 
that, after the integration over arg (x), the linear terms in the boson operators present null coefficients. If we neglect cubic terms 
in the boson operators, we get Hf x written as: 

zjEx zjEx , ttEx 

n l — n Col' n Self' 



H, 



Ex 
Col 



2(2%) 



d 2 q--L l d 2 xV(x)e-'re i ^-^^l 2 



X 2 



Y,V^ G V - ly - l+n (Cf)b n (q) + G V - 1+ny -i(Q*)bl(q) \x{h.c.}, 

n>0 I 



(B9) 



Ex 



1 



Sd f~2%l 2 



_2 /j^2 

d 2 xV(x)e- 1 TLl_ l f - 



n>0 L 



Ly-l 



— Ly-1- 



dqb\ (q) b n (q) . 



Expanding eq.( |B9l ) we will have different exchange effective potentials, v /(q), associated with the different combinations of 
operators b n and b , obtained. For example, the one associated with the combination b n (q) b , (q) is given by: 



0(q) = - 



2 (2n) 2 q 



rf 2 xy(x) e -T g ^-5-)/2 Gv _ 1;V _ 1+;i (Q*)G;_ lv _ 1+n ,(Q*). 



(BIO) 



The integration variable, x, is different from the variable of the integrand, Q*= q* ix*. However, if we decompose the 
function G(Q) using the matrix product (IA2b . apply the Fourier transform, 



e~ % l 2 G m . n (-ix) G* m , n , (-ix*) = — / d 2 ke-^G n , A (k*) G m> (k) e~ k > 2 



2% 

realize the integration over x and, finally, contract the G function using (IA2b again, we get: 

'nn' 1 



(Bll) 



^/^(^/V^-^^Gv^+^^+^k^Gv^^Kk*). 
2(2%) q J 



V 



The angular portion of this integral, k = ke'® k , is given by: 

j rfe k /(c°s%9 y -sin9*fe)g%(>-"') = 2%J n ,_ n (kqje^i ("'-"), 

(B12) 

where we use eq.(3.937) from Ref.[39] and _7„ is the Bessel 
function of order n. Now, due the presence of the coefficient 
e iQ q (n -n) j n eX p ress i on above, only the terms n = n 1 will 
survive after the angular integration over arg (q) = 9 q in ( |B9b . 
Therefore, we will need only 

v^(q) = --f-!^ / dW(k)e-^ 2 (B13) 

k 2 l 2 

XL v _l+„(y)Lv_l(y)^ • 

Repeating the same procedure as above to the other co- 
efficients, we verify that the ones associated with the 



terms b n (q) b n i (q) and b\ (q) b^, (q) will be proportional to 

e ± M n ' +n ), instead of e m "( n '- n ) as in (lBT2l Therefore, once 
n, n' > 0, they will disappear after the angular integration over 
0q in dB9b . The coefficient associated with b\ (q)b n i (q) will 
be the same as the one in dBl3T) . In this way, remembering 
that V(k) = 2ne 2 /ek, we get the expressions ( 13.61 ) to ( 13.9b - 



Appendix C: CALCULATION OF THE FUNCTIONS G< e (r) 

In this Appendix we present the calculation of the Green's 
functions G< e (t), defined in (17.10b . The method used is the 
common one in bosonization procedures. We suggest Ref . [ 1 6] 
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for a review. We have 

w m (e,f)vi4 (0,0) 



= urn — (U m ( t ) e -i^ e -i® m m e ^m(o,o)uV 

e^02%£ \ W m /0 

= Um J_g-'(rtt*+») / e -i© m (e > g i® m (o,o)\ 

e^O 27CE \ / ' 

where we use (17.8b . So, we just have to calculate the propaga- 
tor 

K m e (?) = lim — / e -^(e,0^© m (o.o)\ 

' W e^0 27ie\ /O 

Using the identity = e B e A e c , where C = [A,B] /2 and 
[C,A] = [C,B] = 0, we can write the above expression in a 
normal ordered form, or use directly the relation 

e iih e hh\ = g(^iPi^p2+l(^P?+A|Pi)) , 

/ 

valid for bosonic operators (3 acting on eigenstates of a free 
boson Hamiltonian 18 , and we get: 

1 



where 



K m (Q,t) = lim—e J ^\ 
e^O 27IE 



J m (6,0= D n (9,0 D m (0,0), 



and 



D m (9,0 = (0 m (9,f)0 m (O,O)) 

^0 n=l " 



Here we have used the expression (17.7b of m (9,0- There- 
fore, 

iG>(t) = lim J- e - ! '(^+ 9 ) 
m ' e v ; e^o 27ie 



J5"<**»£(*)E£=i [ e -" l8 ,- /£ "(«)'-i] 



where we use ln(l —x) = — L"=i x " / w > wrtn x = e~ E , in the 
last equation. The calculation of G< e (t) is analogous to the 
one presented above. 



Appendix D: CALCULATION OF THE FUNCTIONS G< e (f ) 

Including the electron spin, the fermionic operator evolves 
in time as 



f m , s (9,0 



1 

hm — U m . s (t)e 

9 -ne/2 



®m,s(Q,t) = i Y — 1=- \ dq<P m (q) 
„ti Vn Jo 

x 'e- in %„Aq,t)-e M btJq,t) . 



Using d5.4| l and d5.5l l we can express the operators b niS (q) 
in terms of the charge and spin bosonic operators, b n (q) p and 
b n (q) a . Then we will have 



®m,s (9,0 



„-ne/2 



— f , 



dq<£>,„ (q) 



x < e 



-m6 



Due the separation of the total Hamiltonian Hj in the charge 
and spin modes (see (15.7b ) we will have the trivial evolution, 

b n (q,t) i = e-^^'bn («?),, with E$ (q) given by £8)- 

Now, the retarded Green's function associated with the in- 
clusion of a particle in the quantum state \n,m,s) is given by: 



iG n%,s (* ) 



Q{t)(c n , m , s (04, m .,(0) 

+ c n,m,i (0) c n,m,s (0 , 



(Dl) 



/■27C 

/9(0 / dQe 
Jo 



inO 



where 



*we(0 = (vL (0,0)^(9,0)- (D2) 



Now, due the commutativity between the operators 

blP (q,t) p and bn (q,t) a and the separation of the total 
Hamiltonian in the form (15.7b . we have the characteristic 
factorization 27 ' 28 of the operators G<^ Q (t) into charge and 
spin propagators. For example: 



iG m , s>e (0 = e-tf'e-BK^pfl (t) K m>afi [t) , 
K m ,i,o(t) = 



1 e ll»=l ^e-^tfd^Me-^W 
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where the factor 1/2 in the exponent comes from the factor 
1 / Vl in the definitions (15.4b and (15.5b of b n (q,t)j ■ Calculat- 
ing iG^ l s % (0 in a analogous way and substituting in (IDlb we 
get (I7.22j ). with 



^m,6 ( f ) = ^m,p,9 ( f ) K m c; Q (t ) 
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